Transport in one dimensional Coulomb gases: Prom ion channels to nanopores 
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We consider a class of systems where, due to the large mismatch of dielectric constants, the 
Coulomb interaction is approximately one-dimensional. Examples include ion channels in lipid 
membranes and water filled nanopores in silicon films. Charge transport across such systems pos- 
sesses the activation behavior associated with the large electrostatic self-energy of a charge placed 
inside the channel. We show here that the activation barrier exhibits non-trivial dependence on the 
salt concentration in the surrounding water solution and on the length and radius of the channel. 



I. INTRODUCTION 

In a number of quasi-one-dimensional systems the in- 
teractions between charged carriers follow (for a certain 
range of distances) the one dimensional Coulomb law: 
$(x) ~ \x\, well known for parallel charged planes. Ex- 
amples include: ion channels in biological lipid mem- 
branes^ 2 ", water-filled nanopores in membranes for de- 
salination devices^, water filled nanopores in silicon oxide 
films^£*&, free-standing silicon nanowires££ and others. 
Their common feature is the presence of a quasi Id chan- 
nel with the dielectric constant greatly exceeding that 
of the surrounding 3d media. As a result, the electric 
field is forced to stay inside the channel, leading to the 
Id interaction potential, mentioned above. The purpose 
of this paper is to study the charge transport through 
such systems as a function of the carrier concentration, 
length, temperature, etc. We show that these systems 
posses very peculiar transport properties, qualitatively 
different from those found in the examples with shorter 
range interactions. 

To be specific, we focus on the water filled channels. 
One example is ion channels in lipid membranes. Such a 
membrane consists of a L = 5 nm thick hydrocarbon layer 
with the dielectric constant K2 — 2, surrounded by water 
with the dielectric constant k\ ~ 80. Due to the large 
ratio K1/K2 — 40 3> 1 the electrostatic self-energy of a 
charged ion inside the hydrocarbon layer is huge, mak- 
ing the pure membrane impermeable for ions from water. 
The only way for ions to cross the membrane is through 
the water filled channels formed by proteins embedded 
into the membrane 1,2 . Radiuses of such cylindrical chan- 
nels, a, vary from 0.3 to 0.8 nm (we are concerned only 
with the passive channels without motors). 

Another example is a water-filled nanopore made in 
a silicon film Such nanopore may have the radius 
a ~ 1 nm and the length L ~ 20 nm. Silicon oxidizes 
around the channel, giving K2 — 4. Thus, also in this case 
k-i/k-2 = 20 3> 1. Yet another example is provided by 
water-filled nanopores, say, with L ~ 40 nm and a ~ 1.5 
nm in cellulose acetate films used for the inverse osmosis 
desalination. In this example K2 — 2 and again K\/k% = 
40 >1. 

Keeping in mind one of these examples, let us con- 



sider a single water filled channel in a macroscopic mem- 
brane or a film separating two reservoirs with salty water 
(Fig. [TJi). We assume that walls of the channel do not 
carry fixed charges. A static voltage V, applied between 
the two reservoirs, drops almost entirely in the channel 
due to the high conductivity of the bulk solution. One 
can measure the ohmic resistance, R, of the channel as a 
function of the temperature, T, concentration of mono- 
valent salt, c (for example KC1) and parameters of the 
channel L, a, k\,K2- The main goal of this paper is to 
evaluate R(T, c, L, a, k±, k%) for the very nontrivial case: 
Ki ^> «2- This is a completely classical {% — 0) problem. 




—I % h 

FIG. 1: Electric field of a charge in a short (a) and a long 
(b) cylindrical channel with the large dielectric constant ki, 
in the membrane with the dielectric constant K2 << Ki. L is 
the channel length, a is its radius. The self-energy barrier is 
shown as a function of coordinate x for the case of the short 
channel. Ul is its maximum height. For a long channel £ is 
the electric field escape length from the channel. 

For a low concentration of salt c the charge transport 
is due to the rare events, when there is a single (e.g. pos- 
itive) ion inside the channel. We assume that the radius 
of the ion, &, is smaller than that of the channel a, so that 
the ion is totally surrounded by water. It is easy to see 
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that the electric field of the ion at distances larger than a 
is deformed because the electric field lines avoid entering 
the media with the small dielectric constant K2 (Fig.Q. 
This leads to the enhancement of the electrostatic self- 
energy of the ion, which creates a barrier U (2) , where 
x £ [-L/2, L/2] is the ion coordinate inside the channel. 
This barrier is the difference between the self-energy of 
the ion at the point x inside the channel and the self- 
energy in the bulk. It is the maximum of the barrier, 
U(0) = Ul, that determines the resistance of the chan- 
nel, Fig. [Hi. 

For a relatively short channel, all field lines are con- 
fined inside the channel until they reach its end and go 
out to the bulk of the water (Fig. QJ,). Therefore, the 
electrostatic self-energy barrier Ul oc L and does not de- 
pend on K2- On the other hand, in a very long channel 
(Fig. nja) the lines eventually escape to the surrounding 
media with the small dielectric constant K2 before they 
reach the exit to the bulk water. As a result, Ul sat- 
urates at some value Uoo, which depends on the ratio 
but does not depend on the length L. 

The self-energy barrier was studied by many au- 
thorsSiifliiiiiiiii^. We shall review their results in sec- 
tion[n] Here we calculate the barrier for the simplest case 
of the short channel^. The electric field Eq at a distance 
x > a from a charge located in the middle of the short 
channel is uniform and given by the Gauss theorem 

Eo = ^- (1) 

The energy of such field in the volume of the channel 
determines the maximum height of the barrier (Fig.^i) 
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where the zero argument indicates that the result is valid 
in the limit of vanishingly low salt concentration. For a 
narrow channel the barrier, Ul{0), can be much larger 
than ksT, making the resistance R exponentially large: 



R ~ exp 



Ul 
k H T 



(3) 



In this paper we study what happens when the con- 
centration of the monovalent salt, c, in the bulk solution 
increases and more ions enter the channel. They first ar- 
rive as compact neutral pairs of oppositely charged ions. 
Indeed, such pairs can enter the channel without the self- 
energy barrier. The one-dimensional character of the 
electric field between the charges creates a strong con- 
finement between them (Fig. Indeed, when a pair of 
oppositely charged ions separates by a distance x, the 
uniform electric field between them, 2Eo , creates the 
confining "string" potential $(x) = e^o|x|. This situa- 
tion reminds two quarks confined in a meson. Condition 
$(x) = ksT defines the characteristic thermal length of 
the pair: 
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ksT kBTnia 2 a 2 
~eE^ ~ 2^2 ~ 21^ 
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FIG. 2: Electric field of a pair of positive and negative charges 
inside the channel. When charges move away from each other 
the energy of the electric field between them grows linearly 
with the separation. 



where Ib = e 2 / (KiksT) is the Bjerrum length (for water 
at the room temperature Ib = 0.7nm). One may think 
about the pair as a one-dimensional object if xt > a or 
a > 21b- We first assume that this condition is fulfilled 
and later return to narrower channels and show that our 
results are valid (with a small modification) for the latter 
case, too. 

At small enough concentration, c, (and/or tempera- 
ture) the thermal pairs do not overlap and the channel 
acts as an insulator with the exponentially large resis- 
tance, cf. Eq.|0). We study how the barrier decreases 
when concentration of pairs grows. At small c we calcu- 
late a relatively small correction to the barrier, which 
however may be larger than ksT and therefore leads 
to the exponential decrease of R with c. In the large 
c limit, one may expect that the thermal pairs overlap, 
screen each other and free individual charged ions. In 
other words, one could expect an insulator-to-metal (or 
deconfinement) transition at some critical concentration. 
This does not happen, however. We show below that, 
due to the peculiar nature of the Id Coulomb potential, 
the barrier proportional to the system's length persists 
to any concentration of the salt, no matter how large it 
is. Its magnitude, though, is exponentially suppressed at 
large salt concentration, leading to the double exponen- 
tial decrease of the resistance with c. The absence of the 
phase transition is in agreement with the earlier studies 
of the thermodynamics of Id Coulomb gasi^iZ*^. 

This paper is organized as follows: Section ITT1 reviews 
the existing literature and formulate our main results for 
the effect of salt on the channel resistance. In Section Hill 
we present qualitative explanations of the results. Sec- 
tions IIVI and [V] are devoted to the quantitative deriva- 
tions for the cases of short and longer channels, corre- 
spondingly. In section IVII we consider modifications of 
the results due to the deviations from the Id model. The 
latter are essential for the very narrow channels. Sec- 
tion EH] treats the dynamics of ions inside the channels, 
that allows to evaluate the pre-exponential factor in the 
resistance. Finally the conclusions and examples are for- 
mulated in section IVlIII 



II. MAIN RESULTS 

Nontrivial physics of the ion transport through the 
narrow channels in the lipid membranes was recognized 
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by ParsegiarA He was the first who studied the barrier 
height L^, in the infinitely long channel. He wrote 



Uoo{0) = ^ F (11 
ma \k 2 



(5) 



and tabulated the function F. Later, numerical solu- 
tions were obtained for a channel of a finite length 101112 . 
Recentljiiii^i^ approximate analytical expressions for 
Ul(0) in the limits of short channel (L <C £) and long 
channel (L 3> £) were derived along with the expression 
for the characteristic length of the channel, £, separating 
these two regimes. 

Let us review the results of Refs. ll3ll4ll1 We already 
presented their result, Eq. for the energy barrier of 
an ion in a short channel [L <C £) . For a very long chan- 
nel (Fig. njj) the electric field of the point charge decays 
exponentially with \x\ due to the leakage of the field lines 
into the media with the small dielectric constant 



E(x) = E Q e 



(6) 



where the characteristic length £ is found as a solution 
of the equation: 



2k 2 J 



(7) 



The exponential decay of Eq. © is replaced by the true 
Coulomb asymptotic E = e/(n 2 x 2 ) : at \x\ ~ £]xi(ki/k2), 
where E -C Eq. 

The self-energy of the point charge in the infinite chan- 
nel may be expressed through £ as: 



£/oo(0) 



e 2 £ 



2 ■ 



(8) 



in agreement with Eq. (J5J). For the protein channels in 
the l ipid m embranes: ki/k 2 — 40 and £ ~ 6.8 a, see 
Refs.l9ll5l 

We augment the literature results for the empty chan- 
nel by showing below that for the channels with interme- 
diate length: 



U L {0) = C/oo(0)tanh 



2£ 



(9) 



This equation correctly matches Eq. © for L <C £ and 
Eq. JSJ for L ^> £. Eq. agrees with numerical calcu- 
lations^ of U L (0) within 10%. 

The goal of this paper is to study effects of a finite 
concentration of a dissociated salt, c, in the bulk solu- 
tion. It is known that for a monovalent salts like KC1 
or NaCl c ~ Co at low total concentration of salt, Co. At 
larger Co > 0.3M (M means a mol per liter) the ratio c/cq 
decreases, but stays close to 0.8 till the saturation limit 
c ~ 6M is reached^. 

First attempts^ to include screening by salt used mean 
field Debye-Hiickel or Poisson-Boltzmann approxima- 
tions. We will show below that discreteness of charge in 



one-dimensional Coulomb case is so important that these 
approximations fail to describe the resistance barrier. 

As explained in the Introduction, the ions easily enter 
the channel in neutral compact pairs. Let us estimate the 
concentration of typical thermal pairs with the arm xt 
inside the channel. First, imagine a channel with k 2 = K>\ 
and the same dimensions L and a as the channel at hand. 
Then the concentrations of ions in the channel is the same 
as outside of it. In such a case the Id concentration of, 
say, positive ions is n = ira 2 c. Using n we can estimate 
the concentration of compact pairs in the channel as: 



2u 2 Xt 



2an , 



where 



nxT 



irca XT 



irca 
~2U 



(10) 



(11) 



is the dimensionless concentration of dissociated salt. In 
Eq. (|10fl at a << 1 the factor 2a may be understood as 
the probability to find a negative ion at the distance xt 
from the positive one. 

In the channel with n\ 3> n 2 we find at small c the 
actual concentration of single ions is much smaller than 
n, so that n is a purely formal construction in this case. 
However, the concentration of pairs n p with the arm less 
than xt is still given by Eq. (fTTTf) . because these pairs 
enter the channel practically without a barrier. As a 
result, the condition of the small overlap of pairs may be 
written as u p xt = 2d 2 <C 1 or just a < 1. Thus, a is 
the main dimensionless parameter of our theory, which 
describes to what extent the pairs overlap. 




a 



FIG. 3: The function / (a) = U L (a)/U L (0) for a short chan- 
nel. Its asymptotic limits valid for small a (Eq. 1121 ) and large 
a (Eq. 1)14^1 are shown by dash-dotted and dashed lines. 

Now we formulate our results starting from a short 
channel, where electric field lines do not escape from the 
channel, L < £. At a <C 1 all ions form compact pairs 
and the effect of interactions between the pairs is small. 
Nevertheless, there is a linear in a <C 1 correction to the 
activation barrier: 



U L {a) - (1 - 4a) U L (0) , 



(12) 



where Ul(0) is given by Eq. J5J. Although the rela- 
tive correction is small, it may have a profound influ- 
ence on the resistance of the channel. Indeed, since 
AaU L (0)/(k B T) = nL, according to Eq. ©: 

R(a < 1) ~ R(0) exp{-ni}, (13) 

where nL = cira 2 L is the number of ions in the volume 
of the channel in the bulk solution (not the actual num- 
ber of ions in the channel). At nL ^> 1, the resistance 
is exponentially decreased in comparison with the zero 
concentration limit. 

In the large c limit, when a > 1, one may expect that 
the overlap of thermal pairs leads to the deconfincmcnt 
transition at some critical concentration, a c — 1. We 
show that such transition does not occur. Instead, the 
barrier proportional to the channel's length exists at any 
concentration of the salt. Its magnitude, however, is ex- 
ponentially suppressed at high salt concentration, i.e. for 
a > 1: 

U L (a) = 2 7 7r- 1 / 2 a 3 / 4 exp{-8 % /^} U L (0) . (14) 

Thus the channel is always an insulator (one may call 
it a swperinsulator, since lni? oc L, rather than R cx L 
as is the case for Ohmic resistors). The a dependence 
in Eq. (|14fl is similar to that obtained in Ref>i& for the 
inverse effective dielectric constant. 

We plot numerically computed function fo(a) = 
Ul(c()/Ul{0) along with its two asymptotic Eqs. iTT2|) 
and (1141) in Fig. [3] It is easy to verify that fo(a) can 
be fitted by fa(a) ~ exp(— 4a) (not shown) within the 
accuracy of 10% in the range < a < 0.94, where /o(a) 
decreases almost 50 times. 

We turn now to the results for longer channels. The 
short channel results, outlined above, are valid for L < 
L c , where L c = L c (a) is a certain threshold length. At 
a <C 1 we show that L c = 7r\/2a£, so the correction 
term in Eq. (|12fl is valid for smaller channel's length. At 
L ^S> L c the barrier saturates at the value 

U oo (a) = (l-Aa 2 ]n(l/a))U oo (0), (15) 

where t/oo(0) is given by Eqs. Q, 

In the «» 1 limit, we obtain the exponentially large 
threshold length L c (a) ~ ^exp{4y / a}. Consequently, 
Eq. I|14|) is only valid as long as L < L c (a). For even 
longer channels, L 3> L c (a), we obtain: 

M«) = 2 5 7r- 5 / 4 a 3 / 8 exp{-4^}[/ co (0). (16) 

Fig. ^] shows the ratio fao(a) = U oo (a)/U oo (0) in the 
infinite channel along with that in the short one. We 
also calculated Jl(o) = Ul(o)/Ul(0) for several different 
ratios of L/£ and plotted these results in Fig.^J 

For the narrow channels, a < Ib, all the results, men- 
tioned above, are still valid, provided that the dimension- 
less concentration a is substituted by the renormalized 
concentration a e //. The latter is defined and evaluated 
in section IVT1 
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FIG. 4: Effect of screening by salt for channels of different 
lengths. Ratios Ul(o) /Ul{0) as functions of a are plotted by 
the dashed line for the infinite channel [f 00(a)) and by the 
dotted line for the short channel (/o(q)). Full lines represent 
the channels with (from top to bottom) L/£ = 4, 3, 2, 1.5, 1. 

Above we have discussed the results pertinent to the 
exponential activation term in the linear resistance. The 
pre-exponential factor is discussed in section IVIII Wc 
have also considered the exponential suppression of the 
differential resistance at large applied voltage, V. For the 
short channel wc obtain the suppression of the activation 
barrier: 

where the applied voltage eV < 2Ul(o). For a long 
channel where L 3> L c (a), Eq. (fTTjl (with Ul — > Uoo) is 
valid as long as eV < e 2 /(K2L). In the opposite limit, 
e 2 I(k.iL) < eV < 2U OCl (a), the long-range nature of the 
3d Coulomb interaction through the dielectric media with 
K2 is important. It leads to the replacement of Eq. i|17|) 
by the Frenkel-Poole law2£ 

mv)~^{ u ^-f B T L ^ n } 08) 

In the next section we present simple qualitative deriva- 
tions of the above results. 



III. QUALITATIVE CONSIDERATION 

A. Concept of the boundary charges 

Consider a relatively short channel where all field lines 
stay inside the channel (Fig.^i). If the average distance 
between the ions is larger than a each ion may be charac- 
terized by a single coordinate x € [-L/2, L/2] along the 
axis of the channel. We deal, thus, with the Id plasma 
interacting through the one-dimensional Coulomb poten- 
tial. Each ion is equivalent to a charged plane perpen- 
dicular to the x-axis with the charge density ±Eq/2tt. 
As discussed in the Introduction, a single negative ion 
in the middle of a channel has the self-energy Ul(0). 



Half of its electric field exits to the right and half to the 
left. One may say that there are image boundary charges 
q = 1/2 and q' = 1/2 on the left and right boundaries of 
the channel correspondingly (hereafter all image charges 
are measured in units of e). One can imagine that these 
charges are provided by the well conducting bulk plasma 
in the reservoirs. 

The concept of the boundary charges, which are not 
supposed to be integer is central for the present work. Let 
us demonstrate the usefulness of the boundary charges 
by a simple example of the ion's energy as function of 
its coordinate x. If the boundary charge at x = —L/2 
is < q < 1, then the boundary charge at x — L/2 is 
q 1 = 1 — q. Electric fields to the left and right of the 
test charge is 2qE and 2(1 — q)Eo, so that the energy 
of electric field U(x;q) = {n 1 /2)E 2 a 2 [q 2 {x + L/2) + (1 - 
q) 2 {L/2 — x)]. Optimizing this energy with respect to q, 
we find q — 1/2 — x/L and for the energy as a function 
x we arrive at 

*<.>-^(M)-H?)><»- 

This gives the parabolic barrier with the maximum value 
£/l(0) in the middle of the channel (see Fig. QJi). The 
maximal barrier corresponds to q = q 1 = 1/2. 

In a similar manner one may consider a pair of a neg- 
ative ion with the coordinate x\ and a positive ion with 
the coordinate x-i inside the channel. They induce two 
boundary charges q and q' = —q. Writing the interaction 
energy of these four charges, one finds that it reaches 
its optimal value, given by U(xi,x 2 ) = 4(|a; 1 2|/£ — 
x 2 2 /L 2 ) Ul(0), at q = X12/L, where xyx = x\ — x%. At 
small separation, \xi%\ <S L the ions attract each other 
with the string potential eEo\x\ — x%\ so that typical dis- 
tance between them is xt- Such a pair can contribute to 
the charge transport only if the positive and negative ions 
separate and move to the opposite ends of the channel. 
Remarkably, the energy barrier they have to overcome is 
given by Ul(0). It is exactly the same as for the single 
ion. The maximum is reached when | cc 12 1 = L/2, while 
q = -q' = 1/2. 

One can repeat this calculation for an arbitrary num- 
ber of pairs inside the channel. One can see that the 
maximum energy state, the system goes through to con- 
tribute to the charge transport, is always reached when 
q = ±1/2 and always has energy Ul(0). We shall refer 
to such a state with the half-integer boundary charge as 
the collective saddle point configuration. Fig.0illustrates 
how the channel with three originally compact pairs ori- 
ented along the external electric field transfers a unit 
charge. It starts from the state with q = —q' w 0, goes 
to the top of the barrier with q = —q 1 = 1/2 and ends 
up again in the state of compact pairs with q = —q 1 « 1. 
The net result of such a process is a transfer of the unit 
charge across the channel. 

From the fact that the energy needed to transfer the 
charge is completely independent of the ion concentration 
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FIG. 5: Three consecutive states of the channel with three 
pairs, (a) Compact pairs with q = —q' ~ 0, (b) the saddle 
point state with q = —q' = 1/2, (c) Compact pairs with 
q = —q' ~ 1. Other saddle points may be obtained e.g. by 
moving the two same sign ions away from each other. 



one may tend to conclude that the same is true regard- 
ing the activation barrier. Such a conclusion is prema- 
ture, however. Indeed, the activation barrier is given by 
the maximum of the free energy and thus includes the 
(negative) entropy contribution. As we explain below, 
the collective saddle point configuration possesses a huge 
degeneracy and thus has a rather large entropy. This 
entropy reduces the activation barrier down to Eq. i|12|) . 

For a more quantitative consideration one needs to dis- 
cuss the exchange of ions between the bulk reservoirs 
and the channel. Although the concentration of single 
ions in the channel is strongly depleted, compact neutral 
pairs enter the channel freely. Thus, typically there are 
many interacting pairs of ions inside the channel. We 
have to discuss separately their equilibrium (minimum 
energy) state and their collective saddle point (barrier) 
state when they contribute to the transport. It is con- 
venient to distinguish the range of relatively low concen- 
trations c, when and the range of high concen- 
trations, a > 1. The concept of the boundary charges, 
introduced above, is helpful in both cases. 



B. Short channel at low concentration of salt 

In this section we consider the case of low concentra- 
tion of ions in the bulk salt solution, such that a -C 1. A 
typical concentration of pairs inside the channel is even 
lower: n p = 2na <C n. To transfer the charge across the 
channel the plasma must go trough the collective saddle 
point state. Such saddle point is characterized by the 
half-integer boundary charge, say q — 1/2; q' = —1/2. 
Half-integer boundary charges apply the external elec- 
tric field Eq across the channel. All dipoles of compact 
pairs turn in the direction of the field, so that positive 
and negative ions alternate. One can see now that in- 
side each pair the electric field changes its sign and turns 
to —Eq (recall the analogy with the uniformly charged 
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planes, producing field Eq in both directions). As a re- 
sult, at the points where charges are located the electric 
field switches between E and —Eq- The electrostatic 
energy of the channel is still equal to E/l(0) and does 
not depend on the positions of the charges. This makes 
the charges free (as long as they are ordered in the alter- 
nating sequence) from being connected into the compact 
pairs. We shall call such peculiar state - the ordered free 
plasma (OFP). 

Since in the OFP state the individual charges are free 
to enter the channel, one expects that the concentration 
of free charges (rather than the much lower concentration 
of the compact pairs) is the same in the channel and in 
the bulk. (Due to the restriction of alternation the actual 
concentration in the channel happens to be half of the 
one in the bulk.) One may say that in the OFP state 
(q = 1/2; q' = —1/2) the gas of the free ions in the bulk 
expands into the channel. This leads to the growth of 
the system's entropy and thus reduction of the activation 
barrier. 

In order to calculate the entropy gain AS of the OFP 
state we assume that from the total number of ions, 2iV, 
some 2AJV enter into the channel's volume AV — ira 2 L. 
The rest 2N — 2 AN ions stay in the bulk volume V. In 
the bulk N — AN positive ions and N — AN negative 
ones produce additive contributions to the entropy. In 
the channel, however, the gas is ordered so that ions of 
different signs alternate and cannot exchange their posi- 
tions. Therefore, the entropy is the same as that of 2A7V 
identical ions. Thus, the total entropy of the system is: 



KB 



AjV)ln 



eV/l 



N - AN 



2 AN In 



eAV 
2A7V/3 ' 



(20) 



where Iq is the normalization volume. The optimal value 
of AN maximizing S is 



A Ar I AV A7 nL 

AN = N= — 

2 V 2 



and thus the entropy gain is: 



AS* = ksnL . 



(21) 



(22) 



Recalling that the free energy at the collective saddle 
point is U L (a) = U L (0) - TAS, we arrive at Eq. 
In this derivation we have neglected entropy change of 
the equilibrium state compared to AS, because n p <C n. 
This simplification is equivalent to omitting corrections 
which are of the order of a 2 in Eq. (|12f> . This calculation 
is justified as long as AN ^> 1. 



C. Short channel at high concentration of salt 

At high concentration of ions in the channel when 
a = nxT > 1 the thermal pairs overlap. Thus, one may 
expect that the ions screen each other and become essen- 
tially free. If this is the case, the concentration of ions 



inside the channel equals to that outside, i.e. the Id con- 
centration is equal to n. It is natural to approach such a 
dense plasma with the mean-field Debye-Hiickel approx- 
imation. Recalling the analogy with the charged planes, 
one needs to calculate the screening of a plane by other 
planes. This is exactly what happens when 3d plasma 
screens a plane. As a result, the electric field of a unit 
charge at a distance x is expected to be: 



E(x) =E e- lxl/xD . 



(23) 



where the Debye length, xd, is given by the standard 
expression: 



xd 



( h B T 
I 87rce 2 



1/2 



k B T 

8ne 2 



1/2 



(24) 



In terms of a one has: Xn ~ n~ x a l l 2 ~ itq" 1 ' 2 , so that 
n^ 1 < xd < xt for a> 1. 

The interaction energy of the system may be written 

as 



W = 



K\a 2 L 



J dE E 2 V{E) 



(25) 



where V(E) is the distribution function of the random 
electric field in the channel. For a ^> 1 it is Gaussian: 

V{E) = (2ir(E 2 )y 1/2 cxp [-E 2 /2(E 2 )] (26) 

with the mean square average value (E 2 ) = AuxdEq — 
W/(ira 2 L). Indeed, fluctuations of the electric field are 
created by the fluctuations of number of positive and 
negative charges in two adjacent segments of the channel 
of the length xd each. A typical fluctuation has opposite 
charges of the order of Q ~ (nxo) 1 ' 2 , which interaction 
energy is fc^T. Thus, it creates the electric field of the 
order of 2QE = 2(nx D )^ 2 E . 

As was discussed above, the activation barrier is asso- 
ciated with the (free) energy difference between the state 
with half-integer and integer boundary charges. Accord- 
ing to the Debye-Hiickel theory any boundary charge is 
completely screened at the Debye length, xr>, an d thus 
is a local perturbation. Then the aforementioned en- 
ergy difference is not an extensive quantity and is pro- 
portional to xd rather than the channel length, L, as 
was the case in the low density limit. Such a qualitative 
change in the scaling of the activation barrier may occur 
only as a result of a phase transition at some a c — 1, 
where the L-dependent part of the activation barrier 
vanishes. It is well-known, however, that the thermo- 
dynamics of Id Coulomb gas, being well described by 
the Debye-Hiickel theory, does not exhibit any phase 
transitionsiSiiLiS. Thus, one either has to explain how 
the Debye-Hiickel thermodynamics may be consistent 
with the L-dependent activation barrier, or find a phase 
transition in the activation barrier which is not seen in 
the equilibrium thermodynamics. We show below how 
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the former of these two alternatives works. There is no 
phase transition in the activation barrier. Instead, the In- 
dependent activation energy decreases exponentially with 
increasing a at a > 1. 

The crucial point is that the Debye-Hiickel approxima- 
tion misses the discreteness of charge. Since the charges 
are discrete, the electric field along the channel may 
change only in the steps of ±2Eo, at the points where 
the charges are located. Therefore at any point in the 
channel, having the boundary charges q and q' = —q, the 
electric field may be equal to 2Eo(q + M). Here M is an 
integer, given by the total charge to the left of the point 
in question. In this manner the information about the 
boundary charges "propagates" through the entire chan- 
nel (in an apparent contradiction to the Debye-Hiickcl 
prediction, Eq. Ipty). 



a AP(E) 



1 


-1 





1 


2 


4. 









FIG. 6: The function AV(E) is the difference of discrete 
distribution functions of electric field for q — 1/2 and q — 
0, which are made of set of ^-functions weighted with the 
envelope function V(E) given in Eq. 1261 . 

Let us show now that this "quantization" of the electric 
field leads to the activation barrier proportional to L, 
though with the exponentially small coefficient. To this 
end we employ the distribution function of the electric 
fields, Eq. (|26[) . as a weighting factor for the quantized 
values of the electric field 2E^(q + M). For example, 
at q = we sum over integers, while at q = 1/2 the 
sum is over half- integers. The transport barrier Ul(o) = 
W(l/2) — W(0) is the difference between these two sums, 
or in other words: 

oo 

U L (a) = ^^ jE?AV(E)dE, (27) 

— oo 

where AT is the difference of discrete distribution func- 
tions for q = 1/2 and q = illustrated in Fig. H3 It 
is clear that Ul{ol) is exponentially small. Performing 
summation over integer M (instead of the integration, 
as in Eq. 1)2 5^) with the help of the Poisson summation 
formula, one finds the activation energy: 

oo 

2 t r 

UL{a) = -^— dEE 2 V(E)cos(7TE/E ). (28) 



As a result, one finds for the activation barrier Ul(o) ~ 
Ul exp{— 7r 2 v / a}. Similarly to the dilute case the ac- 
tivation barrier is proportional to the channel's length. 
Upon increasing the salt concentration the barrier 
monotonously decreases without undergoing the phase 
transition. Comparing the exponent with the exact re- 
sult, Eq. 1)14(1 . one notices that the coefficient is slightly 
overestimated (it 2 instead of 8). The reason is that the 
fields E, contributing to the integral in Eq. J35J), corre- 
spond to the charge of an optimal fluctuation of the order 
Q w enxjj. Therefore both assumptions of the applica- 
bility of the Gaussian distribution and the theory of the 
Debye-Hiickel linear screening to the optimal charge fluc- 
tuations are only marginally correct. The quantitative 
theory of the next section goes beyond these limitations. 
Unlike the activation barrier, the thermodynamics is very 
weakly affected by the "quantization" and is given by the 
Debye-Hiickel theorjiiSiii (up to exponentially small cor- 
rections) . 

D. Long channel 

In the previous sections we dealt only with a short 
channel, where the electric field does not escape from the 
channel to the media with the small dielectric constant 
K2- For finite a the energy of the electric field is de- 
creased due to the screening by other ions. One may 
interpret this decrease of the activation barrier as an 
enhancement of the dielectric constant: k,\ — > K\{a) = 
Ki/fo(a). Indeed, from Eq. @, one may write Ul(cc) = 
e 2 L/(2 Kl (a)a 2 ) = f o (a)U L (0). 

Employing the concept of the concentration- 
dependent dielectric constant, one may formulate 
the results for the self-energy of an ion in the long 
channel. Indeed, both the characteristic length of the 
electric field decay, £ — » and the self-energy 

barrier, C/oo(0) — > Uoo(a), are functions of the dielectric 
constant K\ (see Eqs. (JTJl, 0), which depends now 
on the concentration: n\ — > K\{a). The approximate 
solution of Eq. (7J for is: 

e(a) « a .SSl n V 2 fflM > ) , (29) 

Uoo{a) = — . (30) 

Since according to Eq. fo(a) — e ~ 8 \/" for a > 1, 

one finds «i(a) ~ Kie 8 ^. According to Eqs. 123, CH 
the decay length of the field £(a) ~ e 4v/ " grows exponen- 
tially and the self-energy barrier is E/oo(a) ~ e~ 4v/ ", cf 
Eq. 

For a <C 1 one must take into account deviations of 
E(x) from Eq, cf. Eq. iJBJ. As was discussed in sec- 
tion IIII Bl Eq. I|12|) is based on the fact that the field 
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Eq makes the pairs of positive and negative ions free. 
A somewhat smaller field E(x) leaves most of the pairs 
bound. However, even bound pairs contribute to the en- 
tropy growth. Let us consider a point charge in the mid- 
dle of the channel with the field decaying according to 
Eq. ©• It is easy to see that the average arm of pairs 
grows for the dipoles oriented along the direction of the 
electric field E{x). For the oppositely oriented pairs the 
arm shrinks. As a result, the average pair length be- 
comes 2xt{x) = 2xtEq/{Eq — E 2 {x)] and for the lo- 
cal concentration of pairs instead of Eq. (|10|l we find: 
n p (x) — 2u 2 xt(x) = n p E' 2 l\E' 2 — E 2 (x)\. This estimate 
is valid for \x\ > a£, while at smaller distances n p w n/2. 
Thus the total number of pairs, N p , brought into the 
channel by the central charge may be estimated as: 

L 

AN = 2 J [n p (x) - n p ]dx = 2<aln(l/a), (31) 

<*€ 

The new pairs add the entropic term —UbTAN to the 
free energy of the channel with the fixed central charge. 
This entropic contribution leads to the suppression of the 
transport barrier given by Eq. ijlSfr. 



a long channel one gets similarly to Eq. i|17|) : 

On the other hand, if V/L 3> 2e/{n2L 2 ) the shape of 
the barrier changes drastically. It looses the symmetry 
around x = 0. The saddle point moves to the distance 
I -C L/2 from the left border. The length I is determined 
by the maximum of the sum of the attractive Coulomb 
potential of the charge to the left image and the linear 
potential of the uniform electric field V/L. This maxi- 
mum is located at the distance I — (eL/A^V) 1 ^ 2 . As a 
result, the barrier is reduced by 2eVl/L = (e 3 V/ 'k^L) 1 / 2 , 
thus we arrive at Eq. (|18fl . This result was first derived 
by Frenkcl 20 for the field-induced ionization of Coulomb 
impurities in semiconductors. 

The longer the channel, the smaller is the range of 
applicability of Eq. (|32() and the larger is the range of 
applicability of Eq. l(T%|) . Notice that the applied voltage 
V does not interfere with the effect of the salt on the 
barrier height. In the Eq. JT7|), Eq. and Eq. |JTSJ 
the concentration-dependent height of the barrier deter- 
mines only the voltage at which the activation exponent 
is of order one. Equations 117|) . (|32|) and l|18fl are valid 
only at smaller voltages. 



E. Non-ohmic transport of the channel at large 
voltage 

Until now we dealt with the ohmic transport limited by 
small voltages eV <C ksT. Let us consider the effect of a 
finite voltage on the differential resistance in the opposite 
case, eV > ksT. We start from a short channel. As was 
discussed above, the top of the transport barrier corre- 
sponds to the boundary charges ±1/2. These charges are 
provided by the potential source which makes the work 
eV/2. Thus, the conductance barrier becomes lower by 
this energy and we arrive at Eq. 1)17(1. 

In the long channel the electric field lines of a moving 
charge leave the channel, go to the low dielectric con- 
stant media and eventually return to the water solution 
far from the channel. Thus, each moving charge cre- 
ates the two opposite sign images in the two water-filled 
semi-spaces. These images are similar to the boundary 
charges. Attraction to each of them follows the three 
dimensional Coulomb law. For example, a charge with 
the coordinate x close to —L/2 is attracted to the image 
with the energy -e 2 /[4n 2 (L/2 + x)], if L/2 + x > £. The 
sum of the attraction energies to both water semi-spaces 
has maximum at x = 0. Thus, the top of the barrier 
corresponds to the charge located at the center. In this 
position the image charges are equal q = q' = 1/2. If 
the applied electric field V/L is smaller than the field of 
each image in the center, 2e / (K2L 2 ) , the field only weakly 
affects position of the saddle point. At such conditions 
the voltage only reduces the barrier by eV/2. Thus, for 



IV. THERMODYNAMICS OF THE SHORT 
CHANNEL 

A. Partition function 

In this section we focus on a relatively short channel, 
such that the electric field lines do not escape to the me- 
dia with the small dielectric constant. The calculation 
of the partition function is similar to that in Rcf. 17 . Wc 
briefly reproduce it here to focus on the role of the bound- 
ary charges and to generalize it later to the short channel 
case. The interaction energy of plasma consisting of A*" 
positive and N' negative charges is given by 

1 N+N' 

^=2 E ViVj®(xi ~ Xj) , (33) 
j,i=i 

where Oj = 1 for j = 1, . . . N and ctj = — 1 for j — 
N+l, ... N+N'. The interactions are mediated through 
the Id Coulomb potential: 

<f>(x) = $(0) - eE Q \x\ (34) 

Here $(0) is an artificial constant (twice of the ion's self- 
energy) which is used to impose charge neutrality by tak- 
ing $(0) — » 00 limit. To this end the diagonal terms, 
representing the self-energy of the charges are included 

inEq. iJSJ. 

Following discussion of Section IIII Al one also intro- 
duces two boundary (in general non-integer) charges q 
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and q' at the end points of the channel x — ^fL/2. With 
these two charges included the total energy takes the 
form: 

L/2 

U(q, q') = \ jj dxdx' p(x)<P(x - x')p(x') , (35) 

-L/2 

where the charge density is defined as 

N+N' 

p(x) = cr j 5(x-x j )+qS(x+L/2)+q'S(x-L/2) . (36) 
j'=i 

We are interested in the grand-canonical partition 
function of the gas defined as 



Z L (q,q')= £ 



ol i,(N+N')/(k B T)_ 



1 



N,N'=0 

1 ax 



N\N>\ 



77 (37) 



>< n 



I 2 
'o 



lo 



where p, is the chemical potential of the charges (the same 
for pluses and minuses) and Iq is a microscopic scale re- 
lated to the bulk salt concentration as c — e^/ ksT /Iq. 
Factor na 2 /I 2 originates from the integrations over trans- 
verse coordinates. Since we have in mind the viscous 
dynamics of the ions, rather than the inertial one, the 
kinetic energy is not included in the partition function. 



Formally any configuration of N pluses and N' minuses 
is allowed in this expression. One expects, though, that 
only neutral configurations with N — N' + q + q' = 
provide non-zero contributions due to the presence of 
the large self-energy, $(0). As a result, one expects that 
ZL(q,q') <~ S q ' im - q , where m = N' — TV is an integer 
number. We shall see later that this is indeed the case. 



To proceed with the evaluation of Z^{q, q') we intro- 
duce the resolution of unity written in the following way: 



-L/2 



r ( N+N' \ 

1 = Vp(x) S \p(x) - £ a j 5(x-x j )-qd{x + L/2)-q'S{x- L/2)\ 



J j Vp{x)V6{x) exp 



L/2 
/ L/2 



) (38) 



N+N' 



J dx6(x)p(x)- J2 ^j6{ Xj )-qe{-L/2)-q'e{L/2) 

V-L/2 



Substituting this identity into the expression for the partition function, one finds: 

L/2 



Z(q,q') 



II 



4vr 2 



-oo 
oo 



' Bf JJvp{x)V6{x) cxp< —^-f J J dxdx'p{x)<S>(x-x')p{x')-i j dx9(x)p(x) 

I -L/2 o 



^ Nl 



,nN/(k B T) 



N=0 



2 r \ N co -. / 2 r \ N' 

' N'=0 



I 3 



(39) 



oo f L/2 L/2 

= // e iq0i+i9 ' $f Jv6{x) cxp< J J dxdx'6{x)§-\x-x')0{x') + 2-Ka 2 c J dxco S e(x) 

-oo I -L/2 -L/2 



The integral over e(x) runs over all functions with the 
boundary conditions 9{—L/2) = 6i and 6{L/2) = Of. It 
is straightforward to verify that operator is given 
by <p-\x - x 1 ) = -(2eE )- 1 S(x - x')d 2 x . As a result, 
the functional integral on the r.h.s. of the last expression 



-L/2 

I 

takes the form 

G L {9i,e f ) = IVO(x) e 



(40) 
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where xt and a are given by Eqs. and (|llfl . 

Expression l|40|l represents the "quantum mechani- 
cal" probability to propagate from the point 9i to the 
point Of during the (imaginary) "time" L. The corre- 
sponding "Schrodinger equation" for the "wave-function" 
v & = ^(9, x) has the form: 

Id 2 * x T d^ fA1 . 

-2^- acos ^-"Tfe- (41) 

One may look for the eigenfunctions of this equation: 
^k{9,x) — ^k{9)exp{—2ekx/xT}, labelled by a quan- 
tum number k. The corresponding "Schrodinger equa- 
tion" for the stationary eigenfunctions has the form: 



where m = q — k must be an integer, while \k\ < 1/2; 
i.e. k = [q], where the square brackets denote the frac- 
tional part. We found, thus, that the propagator con- 
serves momentum up to an integer m (umklapp pro- 
cesses): q = q + m. This is a manifestation of the charge 
neutrality. The integer parts of q and q' = —q may be 
arbitrary. Since the integer parts may be easily screened 
by attracting the corresponding number of counter-ions 
from the solution, they do not play an important role in 
the subsequent discussion. One finds then for the Green 
function: 



G x (q,q)=8([q] 



7 (j) 



' x/xq 



(47) 



where we have used the fact that the band dispersion 

acos8 ^k^d) = tk^k{6) , (42) e k ^ s a periodic function of the quasimomentum. The 

common fractional part of the boundary charges plays 
the role of the Bloch quasimomentum. It is important to 
notice that if q is fixed, only the single state from each 
Bloch band contributes to the partition function. 

For a sufficiently long channel, only the lowest Bloch 
band, j — 0, should be retained in this expression, while 
the upper bands provide exponentially small corrections. 
Since the free energy is given by the logarithm of the 



l d 2 ^k(0) 
'2 56» 2 



where tk is the energy. In terms of the stationary eigen- 
functions the propagator takes the form: 



G L (di,6 f ) = **(*<)**(*/) e~ 2 ^ L ' XT . 



(43) 



Finally, the partition function, Eq. I|39l) . is nothing but 
the propagator in the momentum representation and thus 
may be written as: 



Z L (q,q') = G L (q, -q') = £ * fc (g)* fc (-q0 



7 -2e k L/x T 



(44) 

where ^k(q) = f d9/(2n)^ k (9)exp{i9q}. Relation be- 
tween the thermodynamics of the Id Coulomb gas and 
the Mathieu equation (|42|l was realized long ago by 
LenarcU^. He was interested, however, only in the equi- 
librium thermodynamics, rather than the transport char- 
acteristics. 



B. Green function 

One can employ now the well-known results for the 
Schrodinger equation in the periodic potential, Eq. <|42[) . 
to find the propagator, G x (q, q). The wave- functions are 
given by the Bloch-waves: 



propagator, one may also disregard i4°m factors. In- 
deed, the latter contribute to the free energy only a L- 
independent (boundary) term, which is relatively small 
if L xt- Therefore, the propagator may be approxi- 
mated by: 



G x (q,q)^S([q] 



-2e^x/x T 



(48) 



where the (a) is the dispersion relation of the lowest 
Bloch-Mathieu band. Finally, the extensive part of the 
free energy of the channel according to Eq. (|44|l is given 
by: 

or 

F L (q) = -k B T\nZ L {q,q') = fc B T — e<°> , (49) 

XT 

provided [q'] = — [q\. The fact that the free energy is a pe- 
riodic function of q with the unit period reflects the per- 
fect screening of the integer part of the boundary charge. 



Lowest Bloch band 



*f{6) = e 



-ik6„,U) 



(9) = e 



i kit 



E 



U k,m e 



(45) 



where |fc| < 1/2 is the quasimomentum and u^\9) is the 
period-27r (Bloch) function in the j-th band. In the last 
expression one has expanded the Bloch function in the 
Fourier series. The Fourier transform of the eigenfunc- 
tion with the quantum numbers k, (j) is therefore given 
by 



,(i) 

l k,m 



S(q — k — m) 



(46) 



Since the free energy is given by the dispersion relation 
eg (a), let us discuss this function in more details. We 
concentrate on the first Brillouin zone: |q| < 1/2. 

In the extreme low concentration limit, a — > 0, the 
dispersion relation is given by ef^ = q 2 /2. In the dilute 
limit, a < 1, the cosine potential in Eq. I)42|) may be 
treated in the perturbation theory. Save for the points 
q = ±1/2, there are no first order corrections. The sec- 
ond order perturbation theory leads to: 



AO) 



1-Aq 2 



(50) 
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The divergency at q = ±1/2 signifies that the narrow 
interval around these two points, \q ± 1/2| < a <C 1, re- 
quires a special treatment. Indeed, in the narrow region 
around q — ±1/2 points the energy gap opens up. Lin- 
earizing the bare spectrum q 2 /2 in the vicinity of these 
points, one reduces the Schrodinger equation to the spec- 
tral problem for the 2x2 matrix Hamiltonian: 



Sq a 

< I- 



Sq 



(51) 



where Sq = 1/2 ± q. Finding the spectrum, one obtains 
that the lowest Bloch band at \Sq\ < a takes the form: 



c(0) 



1-4^/a 2 



Sq- 



(52) 



,(o) 



At q — ±1/2 the energy reaches its maximum e^, 2 — 
(l-4a)/8. 

In the high concentration limit, a > 1, the amplitude 
of the cosine potential is large. The band is narrow and 
is centered near the ground-state of an isolated (almost) 
parabolic potential well: 



,(o) 



(53) 



The momentum dispersion originates from the tunnelling 
between the adjacent wells. One may treat it in the tight - 
binding approximation and find for the dispersion rela- 
tion: 



-a+^- -2- 4 /o(«)cos(27rg). (54) 



The band- width, 2 3 /o(a), may be calculated in the 
WKB approximation (for the numerical factor see 
Ref.iH): 

2~ 3 /o(a) = 2 4 7r- 1 / 2 a 3 / 4 exp {-8^ } . (55) 



D. Thermodynamics 

One can employ now Eq. I|49|) to extract the thermody- 
namic characteristics of the Id Coulomb plasma (see also 
Ref. IT^) . In particular, we are interested in the salt con- 
centration inside the channel and the average interaction 
energy. The former is given by: 

Up -~2L~dTt~~2Ld^Wt~~ n da ' ( ] 
while the interaction energy is: 



d (F 



----- "sh?;--* 

(57) 

For low concentration, a < 1, and \q =p 1/2 1 > a the 
pair concentration is: n p (q) = 2na/(l — 4q 2 ). This means 



that as long as \q =F 1/2 1 > a, the pair concentration in- 
side the channel is low n p <C n. In terms of n p the free 
energy takes the simple form: F(n p ) — const — Ln p ksT, 
corresponding to the ideal gas of dipole pairs. The in- 
teraction energy is given by W — eE^Lq 2 + ripLksT, 
showing that there is energy ksT per pair (in addition 
to the interaction energy of the boundary charges). 

Consider now the narrow region around q = ±1/2 
points, where \q ± 1/2| < a. Employing Eq. I|52|l . one 
finds for the pair concentration: n p (l/2) = n/2 3> n p (0). 
In fact, it is only half of the naive non-interacting es- 
timate for the pair concentration inside the channel: 
n. This is the OFP state, where the ions are free to 
move as long as they are ordered in the alternating se- 
quence. Indeed, for the free energy one finds F(n p ) = 
const — 2Ln p kBT, indicating an ideal gas of individual 
ions. For q = 1/2 one has W — eE^L{l/2) 2 , meaning 
that the entire interaction energy comes from the bound- 
ary charges and not from the inside pairs. Notice that 
slightly away from the collective saddle point, q = 1/2, 
the plasma in the channel is strongly non-ideal. 

In the dense limit a > 1 the thermodynamic proper- 
ties are almost independent from the boundary charges, 
q. Employing Eq. I|53|) . one finds: n p — n — n/(4 v / a). 
The concentration in the channel is almost equal to the 
bulk concentration. For the interaction energy one finds 
W = nLksT/ (2y/a) in exact agreement with the Debye- 
Hiickel calculation, cf. Eq 1)25(1. Since the interaction 
energy per ion is much less than the temperature, the 
thermodynamics of the plasma is almost identical to that 
of the ideal gas of free ionsi&. Despite of this, the plasma 
possesses a thermodynamically extensive (oc L) activa- 
tion barrier for the charge transport. 



E. Activation barrier 

To calculate the activation barrier, one needs to find 
the free energy of the collective saddle point state. To this 
end we introduce an additional unit charge (say negative) 
at some point, xo, inside the channel. Due to the perfect 
charge neutrality the system must respond by bringing 
two screening boundary charges: q and q 1 = 1 — q, where 
q depends on xq. An example of such dependence is 
discussed above, see Eq. (tT^|l . One can invert the logic 
and claim that by changing the boundary charge q one 
changes the equilibrium position of the test charge, xq. It 
is clear that by varying q from zero to one, the negative 
charge is moved from xq — L/2 to xo = —L/2. The 
activation barrier for such process is determined by the 
maximum of the free energy as function of q. According 
to Eqs. I@3J, (|4*%|l and (|4*§|l the partition function is given 
by: 



Zh{q, =G L (q,q-l) =e 



(58) 



Thus the free energy reaches the maximum at q = 1/2, 
which corresponds to the test charge being in the middle 
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Eq. H62fl ■ It leads to the linear in a correction to the acti- 
vation barrier, which is bound to be smaller than IcbT if 
a < V x t/€- By this reason we shall not discuss it here. 
To proceed we rewrite the Schrodinger equation l|62|) in 
U L (a) =F L (1/2)- F L (0) = U L (0)f o (a), (59) the basis of the Bloch waves, ^^(fl), Eq. gHl- This basis 

diagonalizes the first two terms on the l.h.s. of Eq. I|62|l . 
The only non-diagonal term in this basis is the confine- 
ment potential proportional to 9 2 . In the Bloch basis 
it takes the form -6 k , k > [S j d\ + 2v{ f (k)d k + v^' (k)], 
where v j / (k) = Em"fei 9 fe u fe,m and * = 1,2. Fortu- 
nately, one does not need to evaluate functions v J s ' J (fc). 
Indeed, due to the small coefficient Xy/(8£ 2 ) < 1 in 
front of the 9 2 term, only the highest derivative, 
contributes to the leading order in £,/xt of the WKB ac 
tion. By the same reason one may work only with the 
lowest band, j — 0, and disregard differences between k 
(Bloch quasimomentum) and q (momentum). As a re- 
sult, Eq. (|62[1 takes the form: 



of the channel. The minima are reached at q = and q = 
1, corresponding to the test charge being at x = ±L/2. 
Finally, the activation barrier is given by: 



where the bare barrier Ul(0) is given by Eq. (J2J) and 
the universal function fo(a) is proportional to the Bloch 
band- width: 

/„(a)=8( e f»H- e < 0) (a)) . (60) 

This function is plotted in Fig. 01 The two limiting cases 
of low and high concentration, Eqs. (|12|l and JHJ, follow 
immediately from Eqs. I|52|) and correspondingly. 



V. LONG AND INTERMEDIATE CHANNELS 

A. Partition function 

In the long channel L > £ > xt, where the character- 
istic length £ is given by Eq. (J7J), the electric field lines 
escape to the media with the small dielectric constant, 
k 2 . As a result, the interaction potential gradually looses 
it Id character and crosses over to the 3d Coulomb law 
$(x) ~ 1/(k2X). According to Eq. JfjJ) in the large range 
of distances it has the form: 



$(x) = eE £ e 



-kl/e 



(61) 



Our previous consideration, c.f. Eq. 1|34|) . corresponds to 
the limit £ -> oo. Since $(0) = eE £ = 2U OD (0) (c.f. 
Eq. ||SJ) is now finite, there is no reason to expect the 
strict charge neutrality. It is straightforward to invert 
Eq. J5TJ: $- x (a: - x 1 ) = (2eE )~ 1 S(x - x')(^~ 2 - d 2 x ). 
Thus, one can repeat the mapping onto the quantum me- 
chanics and arrive at the following Schrodinger equation 
for the wave-function ^{9, x): 



1<9 2 * 



2 d9 2 



— a cos 6 ^ 



l -( X A e 2 * = - x ^. (62) 

2 V 2£ J 2 dx v ' 



As before, the grand-canonical partition function, 
Z^{q, q'), of the plasma with the fixed boundary charges 
is given by Eq. ifH|) : 

Z L (q,q') = G L (q,-q>) = ]T * fc (g)* fe (-g') e~ 2 ^ XT . 

k 

(63) 

This is nothing but the Green function of Eq. I|62[l in the 
momentum representation. 

We shall show now that one may use the small param- 
eter xt/(2£) C 1 to find the Green function G x (q,q) 
in the quasiclassical approximation. The results ob- 
tained this way are restricted to the concentration range 
a > v 1 x t /£• F° r l° wer concentrations one may use 
the standard first order perturbation theory in a for 



XT 

2^ 



<9 2 * 



c q 



x T 

'^2~dx 



(64) 



Here ^ = fy(q,x) satisfies the periodic boundary con- 
ditions imposed at q — ±1/2. This is the Schrodinger 
equation for a particle with the heavy mass (2£/xt) 2 ^S> 1 

in the potential e q °' (a). The latter is the dispersion re- 
lation of the lowest Bloch band of the Mathieu equation 
(|42|l . Thanks to the large mass one can employ the WKB 
approximation to evaluate the Green function, Eq. 16311 . 
In this approximation the Green function is given by the 
exponentiated classical action corresponding to Eq. 164|l . 
Such an action may be written in terms of the Lagrangian 
classical "coordinate" q(x) as: 



L/2 

SL{q,q) = — dx 
XT J 

-L/2 



rc2 



t(V x q(x)) 2 +e^ x) 



(65) 



where q(—L/2) = q and q(L/2) = q. A useful analogy 
is that of a Id string with the displacement q(x) and 

rigidity £ 2 , placed in the periodic potential profile e q °\ 
The optimal configuration of such string is given by the 
minimum of the Lagrangian action l|65(l : 



de 



(0) 



= 0. 



(66) 



which must be solved with the specified boundary condi- 
tions. The action is then given by: 



SL(q,q) 



2 

XT 



L£ + t Jp(q)dq 



(67) 



where the canonical classical momentum is given by 
p(q) = y2(e q oS> — £ ) and £ is the integral of motion, 
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found from Eq. 1(66(1 and the boundary conditions. Fi- 
nally, the free energy of the Id plasma with the boundary 
charges q and q' is given by: 

F L (q,q') = -k B ThxG L (q,-q') = k B TS L (q,-q') . (68) 

In the equilibrium state q = q' = 0. The string lies 
in the minimum of the potential: q(x) = and thus 

£ = 6q , while p(q) = 0. The corresponding action is 

Sl(0, 0) = 2Le^/xT and the equilibrium free energy is 
given by: 

F L (0) = k B T — 4 0) . (69) 

Xt 

This is the same result we found for the short channel, 
Eq. (|49|l . Thus the equilibrium thermodynamics is not 
affected by the escape of the electric field from the chan- 
nel (for not too low concentration: a > \Jxt /£)■ On the 
contrary, the activation barrier is essentially different. 



B. Activation barrier 

To calculate the activation barrier consider again a 
negative unit test charge at some point xq inside the 
channel. It induces the screening charges q and q' at 
the boundaries. Since the saddle point is reached when 
the test charge is placed in the middle of the channel, 
we may put q' = 1 — q. Unlike the short channel case, 
there is no strict charge neutrality. Therefore we should 
not expect q to be 1/2 at the saddle point. Instead one 
should consider a function q(x) with the boundary condi- 
tions q(—L/2) = q and q(L/2) = 1 — q, which is a solution 
of Eq. (J2EJ). From the symmetry one expects q(0) = 1/2, 
corresponding to the test charge being at xq = 0. For 
the free energy of the collective saddle point one thus 
obtains: 

F s L p (q) = k B TS L (q,l-q), (70) 

where Sl(9, 1 — q) is the action (|65JI on the optimal string 
configuration and the initial displacement, q, is deter- 
mined from the condition that the string "slides" from 
q(0) = 1/2 down to q(—L/2) = q according to Eq. (iGGl) . 

In the short channel limit, L <C £, the string does 
not bend and thus q(x) — 1/2 at any point. As a re- 
sult £ = and the classical action is given by the 
first term on the r.h.s. of Eq. 1(6 7(1: <Sl(<Z, 1 — q) — * 
5l(1/2, 1/2) = 2Lefj 2 /xT- The activation barrier is thus 
given by Eqs. (|59"|l and JHDJl, as expected. 

In the long channel limit, L>(, the string goes from 
the minimum at q — to the neighboring minimum at 
q = 1, forming the soliton centered at x — Xq. There- 
fore the integral of motion is given by the potential at 
the minima: £ = and the first term on the r.h.s. 
of Eq. (|67|l is nothing but the equilibrium free energy, 



Eq. (|69|l . The activation barrier is then given by the sec- 
ond term on the r.h.s. of Eq. 1)6 7|) . that is the action of 
the soliton: 



Uooia) = U oo (0)8j J2(e i ° ) (a)-e^(a)) dq , (71) 


where 1/^(0) = eE^/2. In the a -> limit tq ] -> <7 2 / 2 
and Eq. (|71|l is trivially satisfied. In general, by analogy 
with the short channel case, one may write U 00 {a) — 

C/oc(0)/oc(a), where f^a) = 8 ^ 2(ef ] - e^) dq is 
the universal function of the dimensionless concentration. 
This function is plotted in Fig. 21 

In the low concentration limit a <ti 1, employing 
Eqs. (|50[1 and JHJ, one arrives at Eq. 1)15(1 . As was men- 
tioned above, this expression is only valid for 1 > a > 
V*r/?- Thus, there is no true logarithmic singularity in 
a — > limit. 

In the opposite limit, a > 1, the band is harmonic 
and exponentially narrow: — = 2~ 3 f a (a) sin 2 nq, 
see Eqs. ((551) . Substituting it into Eq. (O, one 

arrives at Eq. 1(16(1 . Comparing it with the short chan- 
nel result, Eq. 1(14(1 . one finds that the long channel limit 
is applicable when L > £e 4v/ " ^ £. This reflects the 
fact that, according to Eq. 1(65(1 . the characteristic size of 
the soliton grows as an inverse square root of the poten- 
tial height. It is worth mentioning that the applicability 
criterion of the WKB approximation coincides with the 
condition Ul(oi) > k B T. Therefore the results, obtained 
above, are applicable as long as the activation barrier is 
essential. 



C. Intermediate length channel 

We shall discuss now the channel length's dependence 
of the activation barrier. In the limit of zero concen- 
tration, a — > 0, the classical potential takes the form 
of the periodically continued parabola: = q 2 /2 for 
|(7| < 1/2. The corresponding classical action ((67(1 may be 
easily found: S^q, 1 — q) = (£,/2xt) tanh(L/2£), where 
the optimal boundary charge is q = 1/(2 cosh(L/2£)). 
Employing Eq. ijTUj) . one finds Eq. JSJ for the length- 
dependence of the activation barrier. It may be eas- 
ily derived considering the interaction energy, Eq. 133|) . 
with the interaction potential, Eq. ((61(1 . of three charges: 
q and 1 — q at the two boundaries and —1 in the mid- 
dle. Optimizing the energy with respect to q, one obtains 
Eq. ©. 

Yet another derivation of Eq. © may start from con- 
sidering two membrane - water interfaces at x = ±L/2 
as metallic planes creating the series of negative images 
at = (2k — 1)L and positive ones at = 2kL, where 
k is a positive integer. All the images are located along 
the channel axis. The self-energy of the central charge is 
Ul = e0(O)/2, where (f)(0) is the potential of all charges 
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at x = 0. It consists of the potential of the central charge, 
$(0), and the sum of potentials off all the images, <fo(0): 



oo 



<k(0) = 2£; £(0)^(-l) fc e 

k=l 



oL/i + i ' 



(72) 



where we have used Eq. (|61f) . Adding up </>i(0) and 
$(0) and taking into account that £/oo(0) = 2eE^ (c.f. 
Eq. |JSJ), one arrives at Eq. J^J. 

For not very long channel the string does not go 
far away from the maximum of the potential: e 



(0) 
1/2- 



Thus one may expand the potential to the second or- 
der around the maximum: e q 6y 2 — uj 2 8q 2 /2, where 

— 1/2. Substituting this 
form of the potential into the definition of the canonical 
momentum, p(q), and then in Eq. (|67l) . one finds for the 
saddle point action: 



-d 2 e i q ) \ q=1/2 and Sq 



2f (0) T 



2S£ 



(73) 



58 



"1/2 



\J28£ /uj b and 
£ > may be found using the saddle point 



where the integration ran over \5q\ 



approximation. As a result, for L < L c (a), where 



-1/2 



(74) 



the optimal value is 8£ — and the saddle point action is 
S s £ = 2c^) 2 L/xt- This is exactly the same result as with- 
out electric field escape from the channel, cf. Eq. (|5%|) . 
Therefore, deviations of the potential from the ideal Id 
Coulomb law are relevant only for L > L c . 

Employing Eqs. I|52|l and Q54)). one finds for the thresh- 
old length L c (a) = 7r\/2a^ for a <C 1 and L c (a) — 
2£/V7(a) ~ Ce 4x/ " for a > 1. In agreement with Eq. 10, 
L c (0) = and the crossover from the short to the long 
channel is smooth at a — > 0. On the contrary, at a fi- 
nite concentration of salt the length dependence of the 
activation barrier exhibits non-analyticity at L = L c . 
For L < L c the activation barrier is strictly linear in 
L, Eq. I|59(l . while it smoothly saturates to Eq. Ij71|l at 
L > L c . The observables, such as the resistance, do 
not exhibit non-analyticity due to the pre-exponential 
factors. The very similar phenomena takes place in the 
crossover from the thermally activated escape from a po- 
tential well to the quantum tunnelling 22 . 



VI. NARROW CHANNEL 

The previous consideration was based on the model of 
a pure Id motion. The latter is justified if the characteris- 
tic size of the thermal pairs is significantly larger than the 
channel radius: xt S> a. Recalling that xt — a, 2 /{21b), 



one finds that the applicability condition of the Id treat- 
ment is: Ib <C a <C L. These strict inequalities are rarely 
satisfied for real ion channels. The purpose of this section 
is to analyze how the opposite case a < Ib may be incor- 
porated into the theory presented above. The latter con- 
dition may be formulated as fc^T < 6 2 /(ki&), meaning 
that the interaction energy on the scale a is statistically 
significant. 

To take into account the transverse degrees of freedom, 
one may repeat the calculations of section llVl introducing 
the 3d density, p(r), and conjugated field, 9{r). The 
action for the latter takes the form (c.f. Eq. (|40ll ): 



S[9} = 



<f i 



1 



8ttZ, 



■(V ? ( 



2c cos 9 



(75) 



where the 3d integral runs over the volume of the chan- 
nel. The potential field, 9(r), may be parameterized as: 
9(r) = 9q(x) +59(x, g), where g stands for the two trans- 
verse coordinates. The non-uniform part, 59(x,g), may 
be expanded in the basis of the transverse modes of the 
Laplace operator: V^p s = — A 2 ^ s , where <p s (g) have zero 
normal derivative at the channel boundary, and we use 



normalization JcPgtp 2 



7Td 



2 . As a result: 



(x, g) = 9 (x) + ^6> s (x) ip s (g) . 



(76) 



In terms of the new set of fields, 9 S (x) , the gradient term 
of the action (|75|) takes the form: 



L/2 



dx 



-L/2 



(d x 9 o ) 2 +J2{(d x 0s) 2 + X 2 9 2 s 



It is evident from this expression that (9 2 ) = (xtK)- 1 ~ 
Ib/o-. Therefore in the limit a Ib, the transverse fluc- 
tuations may be neglected, 86 -C 1, and the theory re- 
duces to that of the single field, 9o(x), m agreement with 
the previous sections. In the opposite limit, a < Ib, one 
may average the cosine term on the r.h.s. of Eq. I|75[) over 
the fluctuations of 9 S (x) which are governed by the action 
(|77|l . As a result, one obtains the effective action of the 
field 9o(x) of the form of Eq. (|40|) with the renormalized 
amplitude of the cosine term: 

a— b 

a eff = cx T Jd 2 g(e m ) es = 2wcx T f dgge~ u ^^/^ , 



. ( 78) 

where the upper limit of the radial integration is limited 
by the radius of the ion, b. The effective potential in the 
transversal direction is given by 



Ueff(g) 



x lim 



^ A., 



a 

2\x~\ 



(79) 



We have subtracted the self-energy of the ion in the bulk 
solution: e 2 /(2Ki|x|) to regularize the expression at small 



15 



distances, |x| -C a. Indeed, such a local contribution to 
the self-energy is not affected by moving the ion from the 
bulk into the channel. Thus it must be subtracted from 
the chemical potential, fi, to have the proper definition 
of the bulk concentration, c. 

For the cylindrical channel, the effective potential may 
be expressed in the closed form through the modified 
Bessel functions 2 ^: 



oo 

Ueff(Q) = / — 

Kid J 7T 



2 



OO 

E 



Km 



(80) 

The integral on the r.h.s. of this expression is plotted in 
Fig. versus g/a. It diverges as a/[4(a — g)] near the 
channel's boundary, indicating repulsion of the ion from 
the charge image in the media with smaller dielectric 
constant. This divergence plays a minor role in Eq. (J7SJ, 
however, because a — g > b there. 



Ueff 



e^/Kja 




p/a 



top of the entropy gain, discussed in section IIII Bl It 
leads to the increase of the effective concentration a e f / , 
Eq. I|78l) . in comparison with the naive geometrical one, 
a = cxT^a 2 . 

The ratio a e f / jet for three choices of the hydrated 
ion radius b is plotted versus Ib/o, in Fig. [S] One no- 
tices that for a reasonable for NaCl and KC1 valued 
b = 0.281b = 0.2 nm (used in section IVIIII for numerical 
estimates) the effective concentration may be up to 1.6 
times larger than the geometric one (the latter does not 
take into account the finite core radius). The results of 
the previous sections, summarized in section [H] may be 
straightforwardly applied now with the a e / / substituted 
for a. 




FIG. 8: cteff/a as a function of Ib/cl for the three values of 
b/l B . 



FIG. 7: Effective potential in units of e 2 /(«ia) as a function 
of the distance from the axis of the cylindrical channel. 

Remarkably, the effective potential is negative in the 
middle of the channel: U e ft(0) ~ — 1 .le 2 / (nia) . This 
means that the ions have a preference to go inside the 
channel, provided that the long-range part of the self- 
energy, cx L, is not included. As before, the long-range 
part of the interaction energy is described by the field 
9q, which is uniform in the transversal direction. This 
long-range part of the ion's self-energy is related to the 
uniform electric field Eq at distance \x\ ^> a from the 
charge. At smaller distance, \x\ < a the field is almost 
spherically symmetric, the same as in the bulk. As a 



result, the near vicinity of the charge, |x| 



< 



does not 



contribute to the excess self-energy. This "saves" the 
energy l.le 2 /(«ia) per ion. 

Let us return to the collective saddle point with the 
half-integer boundary charges, discussed in section llll Bl 
Such boundary charges produce the electric field Eq, cre- 
ating the OFP state. When an additional ion from the 
bulk enters the OFP state of the channel, its long range 
part of the self-energy is cancelled (by the field Eq of the 
boundary charges), but the energy gain l.le 2 /(Kia) re- 
mains. One may say that each charge reduces the length 
of the string of the field Eq approximately by a, "sav- 
ing" the energy l.le 2 /(«ia). This energy gain comes on 



VII. DYNAMICS 

Having discussed at length the activation barrier and 
its concentration and geometry dependence, we turn now 
to the pre-exponential factor in the expression for the 
channel resistance, Eq. ©. To this end we need a dy- 
namical model for the ion's motion inside the channel. 
It will also allow us to derive non-linear voltage-current 
relations, discussed in section IIII El We shall assume 
that the ions undergo dissipative brownian motion due 
to viscous friction with the water. The latter eventually 
transfers an excess momentum, the ions acquire in the 
external electric field, to the walls of the channel. Under 
these assumptions one may write the Langevin equation 
for the time-dependent coordinate of the i-th ion, Xi{t): 



V dU 

jXi = Oi— - V rji{t) 

L (JXi 



(81) 



where V is an externally applied voltage, 7 is the viscous 
friction coefficient and rji (t) is the Gaussian noise with the 
amplitude given by the fluctuation-dissipation theorem: 



{m(t)T)jtf))=2jk B T6 ij 6(t-t'). 



(82) 
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In the absence of the activation barrier and the interpar- 
ticle interactions these assumptions lead to the channel 
resistance given by the Drude formula: i?o = lL/(2n). 
For a very small concentration a — > 0, the ions cross 
the channel one by one overcoming the activation barrier 
given by Eq. Ijl9(l (but not interacting with other ions). 
The resistance associated with such a process is: 



R — Rq 



H^Z e U L (0)/(k B T) 



(83) 



We shall assume furthermore that the resistance of the 
bulk reservoirs is negligible in comparison with that of the 
channel. This means that the boundary charges adjust to 
an instantaneous configuration of ions inside the channel. 
The boundary charge q(t) is thus found from the mini- 
mum of the interaction energy: dU(q, 1 — q)/dq = 0. We 
have assumed here that there is a single negative excess 
charge, travelling across the channel, and thus q + q' = 1. 
From Eq. I|35[) one then obtains: 



(84) 



This relation imposes the constraint on possible configu- 
rations of ions, Xi. Such a constraint was not included in 
the thermodynamic calculations. To discuss the dynam- 
ics one needs to define the equilibrium constrained free 
energy as: 



(85) 



where the angular brackets denote summation over num- 
ber of ions and integrations over their coordinates with 
the weight given by Eq. I|37|l . Introducing the Lagrange 
multiplier, one obtains the relation between the con- 
strained and the unconstrained free energies: 



iL 

TTXt 



1/2+ioo 



dQ e 



-WrMQ-ifL/m 



(86) 



1/2- 



This integral may be evaluated in the saddle point ap- 
proximation. From the saddle point equation one ob- 
serves that the extensive parts of the constrained and the 
unconstrained free energies coincide at the two extrema, 
q = and q = 1/2. This is the reason why one can dis- 
cuss the activation barrier using the unconstrained free 
energy. The two free energies have distinct non-extensive 
parts, however, which must be kept to determine the pre- 
exponential factor. 

We multiply now the constraint, Eq. I)84[l. by the Drude 
resistance, Rq = r fL/(2n), differentiate it over time and 
employ the Langevin equations l|81[) . As a result: 



Roq{t) = V 



1 x ^ dU 

2^^ ai dx~ 



(87) 



where r](t) = (2n) 1 OtJ7t(i). This is the Gaussian 
noise with the correlator: 

(v(tW)) = J2^{Vi(t)Vj(t')) =2k B TRoS(t-t'). 



It is thus the equilibrium Nyquist noise of the Drude 
resistor Rq. 

Since the boundary charge, q(t), is a weighted aver- 
age of many uncorrelated fast functions, Xi(t), one may 
assume that it possesses a slow dynamics. Thus, one 
may average out the fast motion of individual ions, Xi(t), 
over the equilibrium distribution obtained for an instan- 
taneous adiabatic value of q(t). Averaging this way the 
second term on the r.h.s. of Eq. I|87jl. one obtains^ 



Raq — v 5 — 

oq 



v(t). 



(89) 



The corresponding Fokker-Planck equation for the prob- 
ability distribution function P(q, t) is: 

d t P = -d q I- I=^-[P{V-d q Ft)-k B Td q P}. 

Kq 

(90) 

Following the standard routed, one obtains for the linear 
resistance, R = V/I, in the stationary case, when the 
current / is a constant: 



R = Rq dqe 



,F* L (q)/(k B T) 



dq e 



-F* L {q)/(k B T) 



(91) 



The first integral here is dominated by the vicinity of 
the barrier top, q = 1/2, while the second one - by the 
equilibrium state q = 0, 1. As a result, the resistance is 
given by: 



R = R A L {a)e 



U L (a)/{k B T) 



(92) 



where Ul(o) = ^(1/2) — -Fl(O) is the activation barrier 
and factor Al(oc) originates from the fluctuations around 
the extremal values of the integrals in Eqs. I|86|l and l|91|l . 
Calculating the Gaussian integrals, one obtains for the 
pre-exponential factor: 



A L (a) 



TTX T 1 + LUh 



L UjQLdb 



(93) 



^2,(0), 



where = d^eq°^\ q= o and u\ 

In the dilute limit, a <S 1, from Eqs. (|50"|l and iT^) 
one obtains: ujq w 1 — 4a 2 and u>b « (2a) _1//2 and thus 
Al{ol) w ttxt I \Ly2a) . The divergence at small a is ter- 
minated at a ~ xt / L <C 1 (that is nL ~ 1, so there is 
less than one particle in the channel in average). At such 
a small concentration one can not assume that the q(t) 
is a collective slow degree of freedom. Formally, the inte- 
grals in Eq. I|91|l can not be considered as Gaussian and 
the pre- factor saturates at At (0) = y ~kxt / L as required 
by Eq. ®. 
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In the dense limit, a > 1, one finds that F*(q) « 
F(q) and thus ojq — ujb — %yfo{a)/2 <1. As a result, 
A L {a) = 4x T /(TrLfo(a)) = n^ksT /U L {a). The pre- 
factor grows with the concentration. This growth may be 
associated with the slow motion of the collective variable 
q(t) near the top of the activation barrier. This result is 
applicable as long as Ul(ol) / (fcsT) > 1. 

To conclude, we found that the pre-exponential factor 
A(a), cf. Eq. I|92|) . exhibits the non-monotonous depen- 
dence on the salt concentration. It reaches a flat mini- 
mum at a ~ 0.2, where Al(-2) ~ Sxt/L and increases 
both for smaller and larger concentrations. Employing 
Eq. I|90|> . one may go beyond the linear response theory 
and discuss non-linear voltage-current relation. Since 
such a derivation is relatively standard^ we do not pro- 
vide it here. The corresponding results are listed in sec- 
tions n and mm 

VIII. CONCLUSIONS 

In this paper we suggested a theory of the salt-induced 
reduction of the transport activation barrier in the water- 
filled channels in a media with much smaller dielectric 
constant, than that of the water. Our theory is based on 
one-dimensional nature of the electric field in the chan- 
nel. We have derived a universal dependencies of the 
barrier height Ul(cx) as a function of the dimensionless 
parameter a, given by Eq. flj) . They are presented in 
section^ Remarkably, when the channel radius a grows, 
the one-dimensional character of the problem remains 
valid until at some a the barrier Ul{o) becomes smaller 
than ksT. In this sense one may say that transition from 
one-dimensional channel to the bulk is governed solely by 
a (and not any other parameter, e.g. £/a, or io/a, etc.) 

Now we would like to consider a few examples of ap- 
plication of our theory to water-filled pores. The first 
one is a-hemolysin, a passive protein channel in lipid 
membranes. As we mentioned above, for the water- 
filled channels in lipid membranes £ ~ 6.8 a. For 
the a-hemolysin L = 5 nm, a = 0.75 nm and thus 
£/£ sw 1. For this case, according to Eq. (JSJ we ob- 
tain 11^(0) = 6.8k B T r and according to Eq. © U L (0) = 
0.46£/oo(0) = 3.1fcsT r , where T r is the room temper- 
ature. To describe the effect of salt one can use the 
curve for L/£ = 1 in Fig. with a e ff substituted for 
a. The point where Ul{cx)/Ul(0) = 0.5, corresponds to 
a eff = 0.2. For a = 0.75 nm, l B /a = 0.93 Fig. |S| gives 
a = 0.8a e ff — 0.16, corresponding to c = 0.37 M and 
the bulk concentration of salt Co = 0.46 M. Thus, it is 
enough to have 0.5 M of salt to reduce the activation 
exponent by the factor of two. 

The second example is a water-filled nanopore made 
in a silicon film Such nanopores are considered in 
nanoscience for the transport of DNA, in salty water—. 
We are concerned only with the transport of small ions 



in a e.g. KC1 solution. Let us assume that a = 1 nm 
and the length L = 20 nm. Silicon oxidizes around the 
channel, giving n-i ~ 4 and K1/K2 — 20 3> 1. Using 
Eq. 10, wc find £ = 4.75 a. This gives L/£ ~ 4. For this 
case, according to Eq. (JSJ we obtain C/qo(0) = 3.6fcsT r 
and Eq. @ results in U L (0) = 0.96C/ oo (0) = 3Ak B T r . 
Effect of salt is given by the top full curve in Fig. 0] 
corresponding to Z/£ = 4. In this case, Ul{oc)/Ul{Q) = 
0.5 at a e ff = 0.4 or at a = 0.36 (as we see in Fig. |S] 
for Is/a = 0.7 the ratio a e ///a = 1.12). This requires 
c = 0.26 M and c = 0.31 M. 

The third example is a water-filled nanopore with L = 
40 nm and a = 1.5 nm in cellulose acetate film used for 
inverse osmosis desalination. Let us take n-i ~ 2 so that 
i ~ 6.8a ~ 10 nm , L/£ ~ 4, U L {0) = 0.96f7oo(0) ps 
3kBT r . Effect of salt is again given by the top full curve 
in Fig. Hand U L {a)/U L (0) = 0.5 at a eff = 0.4 which 
in this case corresponds to a = 0.38. This happens at 
c = 0.69 M and c = 0.73 M. 

The last example, we want to consider, is the 
gramicidin- A-like channel in a lipid membrane. In this 
case we use L = 2.5 nm, a = 0.3 nm and obtain 
L/£ = 1.2. According to Eq. JSJ the barrier ^(0) = 
17k B T r and according to Eq. @: U L (0) = 0.55 [/^(O) = 
9.4fcsT r . We see that, because of the very small width, 
we deal with much higher barriers. On the other hand, at 
such a small a even the concentration of salt cq close to 
the saturation limit of 6 M and corresponding c = 5.2 M 
we can reach only a = 0.06 or a e // = 0.084 (for this case 
Is/a = 2.3 and according to Fig. |S| a e ff/a = 1.4). At 
this a e ff interpolating between curves for L/£ = 1 and 
L/£ = 1.5 we arrive at U L (a) = 0MU L (0) = 7.8k B T r . 
Thus, the effect of salt on a narrow channel is weak. This 
could be anticipated from Eq. (|11|) . if we rewrite it as 
a = tt(cIb) 3 (cl/Ib) 4: and take into account that cl B < 1, 
even at the saturation limit Co = 6 M. Contrary to the 
case where cl/Ib > 1 and the range a > 1 is accessible, 
in the opposite limit a/Is ^5 1 the parameter a is bound 
to be small2&. 

In this paper we considered only channels without 
charges fixed on their walls. Methods used here can be 
easily applied for channels with fixed charges of one or 
both signs. We will address such channels in the next 
publication. 
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To carry out this averaging one parameterizes the ions co- 
ordinates as Xi(t) = <ji(q(t) — l/2)/(2n) +5xi(t), where the 
constraint 1841 is reduced to aiSxi(t) = 0. The second 
term on the r.h.s. of Eq. 1871 may be written now as: 



dU 
dq 



_dU 



{* 4 } 



dU 

dxi 



dxi 

~8q~ 



dU 



2n 



where we employed that dU/dq\{ x .y = 0. Av- 
eraging over the equilibrium distribution func- 
tion Z" 1 exp{-U({6xi},q(t))/(k B T)}, one finds: 
(2n)- 1 (j: i ^dU/dx l \ q ) = (dU/dq\ {SXt} ) = dFl(q)/dq. 
Another possible application of our theory is the trans- 
port in single-crystal silicon nanowires- suspended in air 
between two conductors. They have a ~ 1 nm and lengths 
up to several microns. In this case, n\ — 11 and Ki = 1 
so that the ratio K1/K2 is relatively large. We will address 
this application in a future publication. 



